Supplementary Material 



The following tables report parameters and compare vibrational band origines as discussed in the paper Algebraic- 
matrix calculation of vibrational levels of triatomic molecules by T. Sedivcova-Uhlfkova, Hewa Y. Abdullah, 
and Nicola Manini. 



TABLE A. 


Fitting parameters (in 


Hartree atomic units 


for the R i p expansions in Morse coordinates, Eq. (11) of the main 


text, for the studied molecules. The fits are carried out 


in the following 


ranges: _Roc = 1-32 


-4-4.7 a.u., Res = 2.4 


-4-3.78 a.u., 


i?HC = 1.5 


-r 3.5 a.u., _Rcn = 1.6-7- 


3.0 a.u., ii N o = 1.9 -7 


3.2 a.u., Roh = 


1.35-7-3.2 a.u. 






system 


bo 


bi 


b 2 


h 


bi 


bs 


^oc 


0.457685 


0.169226 


-0.021905 


0.015526 


-0.013160 


0.015850 


rt oc 


0.209479 


0.154873 


0.008480 


0.007549 


-0.004657 


0.005804 


#cs 


0.338904 


0.111322 


-0.018657 


0.012254 


-0.011645 


0.008461 


^cs 


0.114859 


0.075453 


-0.000423 


0.004164 


-0.003489 


0.002530 


D — 1 


0.496627 


0.254217 


0.003579 


0.013680 


-0.015465 


0.026987 


^HC 


0.246653 


0.252068 


0.067223 


0.020238 


-0.000780 


0.008569 


-^CN 


0.458835 


0.171067 


-0.020751 


0.019016 


-0.016428 


0.007374 


D— 2 
n CN 


0.210545 


0.157132 


0.009775 


0.008909 


-0.005792 


0.002694 


^ON 


0.445632 


0.118759 


-0.025096 


0.007967 


-0.021604 


0.026137 


r>— 2 

"on 


0.198628 


0.105467 


-0.009673 


0.004369 


-0.009938 


0.011674 


^HO 


0.551962 


0.261286 


-0.003594 


0.003427 


-0.027048 


0.041423 


%o 


0.304789 


0.286850 


0.059773 


0.013633 


-0.007170 


0.014774 



TABLE B. The values of potential-energy parameters for HCN (Ri = 2.0135 a.u., ai = 0.9727 a.u. -1 for CH, and 

7?2 = 2.1793 a.u., a 2 = 1.2290 a.u. -1 for CN), in Hartree atomic units. The form of the potential energy function is given by 

Eq. (3) in the main text. 



fcl 


fc 2 


k:t 


ak 1 k 2 k 3 


fci 


k 2 


k 3 


a.k 1 k 2 k 3 


fci 


fc 2 


fc 3 


dk 1 k 2 k 3 








1 


0.059318 


1 


1 





-0.010473 


3 





1 


0.008907 








2 


0.011186 


1 


2 





0.007719 





1 


1 


0.060746 








3 


0.008006 


1 


3 





-0.002693 





1 


2 


0.003288 








4 


-0.008333 


2 


1 





0.004400 





1 


3 


-0.015908 


1 








-0.000455 


2 


2 





-0.010019 





2 


1 


-0.028872 


2 








0.211513 


3 


1 





-0.005224 





2 


2 


0.005894 


3 








0.009211 


1 





1 


0.016978 





3 


1 


0.037804 


4 








0.030108 


1 





2 


0.019794 


1 


1 


1 


0.022503 





1 





-0.000791 





3 


1 


0.002993 


1 


1 


2 


-0.018163 





2 





0.397690 


2 





1 


-0.013440 


1 


2 


1 


-0.007818 





3 





-0.022029 


2 





2 


0.007265 


2 


1 


1 


-0.021432 





4 





0.018340 
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TABLE Ca. Comparison between observed and computed purely vibrational m = band origins [cm 1 ] of OCS. 



Vl V 2 V3) 


E° hs 


E^ 1 


7riobs TTical 


100) 


858.97 


858.94 


0.03 


2 0) 


1047.04 


1047.05 


-0.01 


200) 


1710.98 


1710.85 


0.13 


12 0) 


1892.23 


1892.41 


-0.18 


001) 


2062.20 


2061.98 


0.22 


040) 


2104.83 


2104.69 


0.14 


300) 


2555.99 


2555.73 


0.26 


2 2 0) 


2731.40 


2731.65 


-0.25 


101) 


2918.11 


2917.90 


0.21 


140) 


2937.15 


2937.17 


-0.02 


21) 


3095.55 


3095.27 


0.28 


060) 


3170.64 


3170.38 


0.26 


400) 


3393.97 


3393.57 


0.40 


320) 


3564.48 


3564.73 


-0.25 


240) 


3762.83 


3762.90 


-0.07 


2 01) 


3768.50 


3768.26 


0.24 


121) 


3937.43 


3937.42 


0.01 


160) 


3990.11 


3989.93 


0.18 


002) 


4101.39 


4101.20 


0.19 


041) 


4141.21 


4140.86 


0.35 


500) 


4224.87 


4224.38 


0.50 


080) 


4242.55 


4242.26 


0.29 


420) 


4391.41 


4391.63 


-0.22 


340) 


4585.25 


4585.38 


-0.13 


301) 


4609.85 


4609.45 


0.40 


2 21) 


4773.22 


4773.36 


-0.14 


260) 


4805.16 


4805.03 


0.13 


102) 


4953.88 


4953.75 


0.13 


141) 


4970.43 


4970.36 


0.07 


600) 


5048.62 


5048.13 


0.49 


180) 


5049.14 


5048.83 


0.31 


2 2) 


5120.98 


5120.64 


0.34 


61) 


5196.01 


5195.72 


0.29 


5 2 0) 


5212.14 


5212.31 


-0.17 


010 0) 


5319.02 


5318.84 


0.18 


440) 


5401.48 


5401.58 


-0.10 


4 01) 


5444.96 


5444.45 


0.51 


3 21) 


5602.47 


5602.64 


-0.17 


360) 


5616.07 


5615.95 


0.12 


241) 


5792.00 


5792.15 


-0.15 


202) 


5801.91 


5801.68 


0.22 


2 8 0) 


5852.00 


5851.66 


0.34 


700) 


5865.18 


5864.82 


0.36 


12 2) 


5959.33 


5959.37 


-0.04 


161) 


6012.06 


6011.95 


0.11 


6 2 0) 


6026.59 


6026.74 


-0.15 


1 10 0) 


6112.59 


6112.28 


0.31 


003) 


6117.58 


6117.64 


-0.06 


042) 


6154.70 


6154.36 


0.34 


540) 


6212.12 


6212.16 


-0.04 


081) 


6257.83 


6257.74 


0.09 


501) 


6273.13 


6272.51 


0.62 


012 0) 


6398.77 


6398.83 


-0.06 
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TABLE Ca. Continued 



(Vl V2 V3) 


£iobs 




TTiobs TTical 
£L/ V — Ej v 


(4 6 0) 


6418.68 


6418.46 


0.22 


(4 21) 


6429.10 


6429.31 


-0.21 


(3 41) 


6611.11 


6611.31 


-0.20 


(3 2) 


6640.10 


6639.80 


0.30 


(3 8 0) 


6650.90 


6650.55 


0.35 


(8 0) 


6674.46 


6674.41 


0.04 


(2 2 2) 


6791.54 


6791.81 


-0.27 


(2 61) 


6823.87 


6823.86 


0.01 


(7 2 0) 


6834.69 


6834.89 


-0.20 


(210 0) 


6903.16 


6902.78 


0.38 


(103) 


6966.17 


6966.32 


-0.15 


(142) 


6980.96 


6981.00 


-0.04 


(6 4 0) 


7017.13 


7017.11 


0.02 


(181) 


7060.66 


7060.59 


0.07 


(6 01) 


7094.22 


7093.59 


0.63 


(0 2 3) 


7123.38 


7123.16 


0.22 


(112 0) 


7179.14 


7179.02 


0.12 


(0 6 2) 


7198.75 


7198.61 


0.14 


(5 6 0) 


7220.13 


7219.81 


0.32 


(5 21) 


7245.75 


7246.00 


-0.25 


(0101) 


7324.96 


7325.20 


-0.24 


(441) 


7423.79 


7423.92 


-0.13 


(4 8 0) 


7445.74 


7445.45 


0.29 


(4 2) 


7472.18 


7471.77 


0.41 


(9 0) 


7476.38 


7476.94 


-0.56 


(0 14 0) 


7480.69 


7481.16 


-0.47 


(3 2 2) 


7616.85 


7617.19 


-0.34 


(3 61) 


7631.94 


7632.02 


-0.08 


(8 2 0) 


7636.35 


7636.71 


-0.36 


(310 0) 


7690.33 


7689.95 


0.38 


(2 4 2) 


7797.96 


7798.39 


-0.42 


(2 3) 


7811.95 


7811.94 


0.01 


(740) 


7816.43 


7816.39 


0.04 


(2 81) 


7859.95 


7859.91 


0.04 


(70 1) 


7908.22 


7907.66 


0.56 


(212 0) 


7957.27 


7957.12 


0.15 


(12 3) 


7958.04 


7958.28 


-0.24 


(16 2) 


8011.27 


8011.42 


-0.15 


(6 6 0) 


8015.70 


8015.31 


0.39 


(6 21) 


8057.05 


8057.29 


-0.24 
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TABLE Cb. Comparison between observed and computed m = 1 vibrational band origins [cm x ] of OCS. 



Vl V 2 V3) 


E° hs 




7riobs zrical 


010) 


520.42 


520.70 


-0.27 


110) 


1372.46 


1372.85 


-0.39 


030) 


1573.37 


1573.52 


-0.16 


210) 


2218.03 


2218.45 


-0.42 


130) 


2412.12 


2412.47 


-0.35 


011) 


2575.31 


2575.30 


0.01 


5 0) 


2635.59 


2635.59 


0.00 


310) 


3057.09 


3057.48 


-0.39 


230) 


3245.26 


3245.70 


-0.44 


111) 


3424.14 


3424.33 


-0.19 


15 0) 


3461.58 


3461.71 


-0.13 


31) 


3615.34 


3615.23 


0.11 


070) 


3704.77 


3704.68 


0.09 


410) 


3889.61 


3889.93 


-0.32 


330) 


4072.69 


4073.15 


-0.46 


211) 


4266.32 


4266.61 


-0.29 


2 5 0) 


4282.91 


4283.09 


-0.18 


131) 


4450.76 


4450.96 


-0.20 


170) 


4517.94 


4517.88 


0.06 


012) 


4607.11 


4607.06 


0.05 


51) 


4666.09 


4665.96 


0.14 


510) 


4715.54 


4715.79 


-0.26 


090) 


4779.23 


4779.16 


0.07 


430) 


4894.35 


4894.77 


-0.42 


350) 


5097.05 


5097.26 


-0.22 


311) 


5104.23 


5104.50 


-0.27 


2 31) 


5280.54 


5280.93 


-0.39 


270) 


5327.03 


5326.96 


0.07 


112) 


5452.46 


5452.68 


-0.22 


151) 


5488.78 


5488.90 


-0.12 


610) 


5534.81 


5535.04 


-0.24 


190) 


5579.46 


5579.32 


0.15 


032) 


5634.23 


5634.08 


0.16 


530) 


5710.14 


5710.51 


-0.37 


71) 


5724.78 


5724.76 


0.02 


0110) 


5857.57 


5857.68 


-0.11 


4 5 0) 


5908.98 


5909.11 


-0.14 


411) 


5932.73 


5933.00 


-0.27 


3 31) 


6104.52 


6104.97 


-0.45 


370) 


6131.88 


6131.81 


0.07 


212) 


6290.91 


6291.31 


-0.40 


2 51) 


6307.12 


6307.37 


-0.25 


710) 


6347.35 


6347.66 


-0.31 


290) 


6376.32 


6376.15 


0.17 


132) 


6466.15 


6466.44 


-0.29 


630) 


6520.02 


6520.33 


-0.32 


171) 


6534.40 


6534.51 


-0.11 


013) 


6615.85 


6615.96 


-0.10 


1110) 


6644.69 


6644.72 


-0.02 


5 2) 


6673.76 


6673.71 


0.06 


5 5 0) 


6714.77 


6714.83 


-0.06 


511) 


6755.61 


6755.84 


-0.23 


091) 


6789.57 


6789.84 


-0.27 


431) 


6921.91 


6922.20 


-0.29 


390) 


6933.04 


6933.14 


-0.10 


013 0) 


6938.60 


6939.16 


-0.55 


3 51) 


7116.51 


7116.93 


-0.42 


312) 


7126.90 


7127.26 


-0.37 
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TABLE Cb. Continued 


(vi Do Vi) 


pobs 


peal 

V 




(8 1 01 


7153.07 


7153 60 


—0.52 


(4 701 


7169.47 


7169.43 


0.04 


(2 3 2) 


7292.42 


7293.00 


-0.58 


(7 30) 


7323.89 


7324.20 


-0.31 


(2 7 11 
V^ 1 ^/ 


7340 06 


7340 35 


—0.28 


(2 11 0) 


7429.19 


7429.58 


-0.39 


(1131 


7457.43 


7457.87 


-0.44 


(15 21 


7493.10 


7493.44 


-0.34 


(6 5 0) 


7515.15 


7515.15 


0.01 


(6 1 11 


7572.02 


7572.22 


—0.20 


(1911 


7585.81 


7586.33 


-0.53 


(0 3 3) 


7630.11 


7630.10 


0.01 


(1 130) 


7712.37 


7713.60 


-1.23 


(0 721 


7722.39 


7722.74 


—0.35 


(5 70) 


7725.25 


7725.09 


0.16 


(5 31) 


7737.80 


7738.39 


-0.59 


(0111) 


7858.90 


7860.06 


-1.17 


(4 51) 


7925.38 


7925.82 


-0.44 


(412) 


7951.73 


7952.04 


-0.32 


(910) 


7951.89 


7952.94 


-1.05 


(4 9 0) 


7958.81 


7959.27 


-0.47 
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TABLE D. The observed vibrational band origins [cm x ] of HCN are compared to those computed using the present 
approach and to those computed in the Reference work by van Mourik et al. [26] . A few predicted levels are also listed. 



(Vl V2 V3 ) 




rrical 


/TTobs ijical 


zrrRef 


TT'obs zjiRef 


\\J z UJ 


1/111 A O 

1411. 4z 


1 A 1 O A 

141Z.4 


1 n 
— l.U 


1 a 1 /i n 
1414. y 


cc 
—0.0 


(UU 1J 


onnc otc 
zUyb.oo 


1 ni t 
Z^lUl. / 


a n 

— 4.y 


01 nn c 
ZlUU.O 


O 7 

—6. ( 


(U 4 UJ 


nono rift 

zoUz.yb 


oon7 
zoU / .0 


A Q 
— 4.0 


oom a 
zoUl.4 


1 7 
1. ( 


(1 U UJ 


QQ 1 1 A O 

ooll.4o 


OOl4.Z 


O 1 

— A. ( 


oon7 7 


O Q 
O.O 


(U z 1J 


O K AO 1 O 

ooUz.lz 


ccno 


a 
— D.Z 


c 1 n n 

ooiu.y 


O O 
— O.O 


(U UJ 


A 1 TA d 1 

41 ( 4.bl 


A 1 OQ Q 
4100.0 


Q 7 
— O. ( 


A 1 7fi O 

41 /D.z 


1 £. 
— I.O 


(U zj 


A 1 TO f!7 

41 ( 6.\j( 


A 1 oc n 

4i5o. y 


IOO 

— Iz.o 


/1 1 1 /I 

4101.4 


O O 
— O.O 


(1 z UJ 


AdQ A Q1 

4bo4.ol 


/ifio^ n 
4DoD.y 


a. 

— Z.D 


/IfiOfi 

4DOD.Z 


1 O 

— i.y 


(U 4 1J 


4ooo.UU 


/I Qf\7 O 

4oy / .0 


n o 

— y.o 


/i om 7 


7 
—0. ( 


(1 U 1J 


Ooyo. / U 


k /i n 1 cc 
04U1.0 


7 Q 
— I.O 


0oy4.4 


n 7 
— U. < 


(V o UJ 


OozO.ol 


0040. D 


1 n q 

— iy.o 


OOO ( . ( 


1 1 n 

— 11. y 


(U Z z!J 


r r-71 on 


OOOO. 


— 11.0 


K C 0*S C 

OOOD.O 


1 a a 
— 14. b 


(1 4 UJ 


bUob.yb 


DU4Z.D 


CC iS 

— O.D 


finoo 7 



6.6 


(U U oj 


0000 £?n 


Dz44. / 


— lD.1 


(SO A O /I 

0Z4Z.4 


IOO 

— lo.o 


[V 1 J 


Dz04.oo 


uz ( U.l 


1 K 7 
— 10. / 


DZDU.O 


1 

— 0.1 


(z U UJ 


DOly.Dl 


con 1 
DOzy.l 


n c 

—y.o 


(SCC1 cc 
DOlo.O 


a 1 
b.l 


(1 Z 1J 


O/Ol. OO 


u (Do. / 


7 /I 

— / .4 


D/DO.O 


7 
— ( .z 


(U 1U UJ 


DoOO.OO 


Doo4. / 


on 

— zy.z 


Do ( y.D 


/I 1 
— z4.1 


(U4 2J 


c;nc: 1 (So 


oyoo. / 


— lo.U 


oyou.y 


n 

— y.z^ 


(1 UJ 




7QQ1 1 




70(sn 1 

/ oby.i 




(1 U zj 


7/1 EC /I O 

f 400. 4z 


*7A i^Q O 


— 10.4 


7/1 (S 1 


a 

— b.z 


(U 8 1J 








7/S 1 7 O 

/bl / .Z 




(U z oj 




Ti^oo n 
r Doo.U 




7iS /II O 

( 041.Z 




(z z UJ 


TO CO K 1 
fOOO.Ol 


TOd 1 T 
f OOl. / 



— 5.Z 


70 C K O 
/OOO.O 



— 1.6 


^1 4 1J 


o 1 A7 ri7 
olU / .y / 


Q1 1 Q O 

ollo.o 


1 n o 
— lU.o 


1 1 n 
ollU.z 



— z.z 


(U Iz UJ 




oonn 




1 n7 c 

oiy ( .0 




(U U 4J 




oooiS n 

OZOD.U 




OOOO 

ozoo.o 




(U z J 




OOOO (S 
OOOZ.D 




OOOO cc 
OOZO.O 




(A U 1 J 


0O0O.O0 


oDUU.4 


1 /I Q 

— 14. 


O CCO /I 7 

O0o4. ( 


n 

u.y 


(1 o UJ 




Q7fl1 O 




oeni /s 

ooyi.o 




(1 z zj 


oolb.UU 


ooon n 

oozy.y 


ion 

— lo.y 


ooon 
oooU.z 


1 A O 

— 14. z 


(U 1U 1 J 




qoc7 n 

oyo < .u 




once A A 

oyo4.4 




(U 4 oj 


oyyo.zz 


nni /i a 
y(J14.4 


— iy.io 


nnnn n 

yuuy.u 


1 O Q 

— lo.o 


(z! 4 UJ 




HI l ~7 r 7 A 

yi / (A 




ni ct a n 

yio4.u 




(1 1 J 




a a c 1 
y401.o 




C\A AC\ 1 

y44U.l 




(U 14 UJ 




y4yu. ( 




n/i 00 ec 
y4oo.O 




(1 U oj 


O A ClR A A 

y4yo.44 


n 1 a n 

yoio.u 


1 n 

— iy.o 


nccno n 

youo.y 


— lz.0 


U UJ 


yoz / .uy 


OlS At: A 
yt)40.4 


IOO 

— lo.o 


nfi 1 n 

ybiy.z 


7 

* .y 


^U z 4J 




yo/ 1. ( 




yo ( 4.0 




(U o zj 




yoou.u 




yb ( o.z 




(z z 1 J 


nni /i a 

yyi4.4 


nnoo 

yyzo.z 


IOO 

— lo.o 


nnoo n 

yyzz.y 


cc 
—0.0 


(1 1U UJ 




1 nnni n 
1UUU1.U 




nnno n 

yyyo.y 




(1 4 zj 




1 ni 7^0 7 
1U1 (6.1 




1 m a c *s 
lUlbo.b 




(U Iz 1 J 




1 no^7 n 
lUzD / .U 




1 noiSiS /i 
lUzOu.4 




(U U Oj 




1 noniS n 
lUoUo.y 




1 nQfi/i 1 
!UoU4.1 




(0 6 3) 


_ 


10373.5 


_ 


10364.9 





(2 6 0) 




10476.1 




10460.5 




(2 2) 


10631.4 


10652.2 


-20.8 


10636.7 


-5.3 


(016 0) 




10751.0 




10749.7 




(181) 




10766.2 




10758.0 




(12 3) 




10870.0 




10871.1 




(3 2 0) 




10937.3 




10925.3 




(010 2) 


10974.2 


11006.9 


-32.7 


11006.5 


-32.3 
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TABLE D. Continued 



(«i V'l V3) 
(0 4 4) 
(2 41) 
(112 0) 
(162) 
(0141) 
(104) 
(3 01) 
(0 2 5) 
(0 8 3) 
(2 8 0) 
(2 2 2) 
(018 0) 
(1101) 
(143) 
(3 4 0) 
(0 6) 
(012 2) 
(0 64) 



iobs 



11015.9 



11516.6 
11674.5 



11040.8 
11239.0 
11277.5 
11500.5 
11542.9 
11550.9 
11684.5 
11698.5 
11714.2 
11756.7 
11974.5 
11978.0 
12060.1 
12206.2 
12213.4 
12307.5 
12310.6 
12392.9 



-ical 



-24.9 



-34.3 
-10 



E Roi 
11035.7 
11225.9 
11271.7 
11489.4 
11536.1 
11549.6 
11672.8 
11688.1 
11710.1 
11744.5 
11969.9 
11977.7 
12055.7 
12192.5 
12200.6 
12304.5 
12311.7 
12382.0 



pRcf 



-19.8 



-33.0 
1.7 
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TABLE E. Comparison between observed and computed vibrational band origins [cm *] of NO2. 



Vl V 2 V3) 


E° hs 


E^ 1 


rpobs Trical 


010) 


749.649 


749.48 


0.17 


100) 


1319.794 


1320.67 


-0.88 


2 0) 


1498.34 


1498.22 


0.12 


001) 


1616.852 


1618.34 


-1.49 


110) 


2063.118 


2062.35 


0.78 


030) 


2246.04 


2245.90 


0.14 


011) 


2355.151 


2355.65 


-0.50 


200) 


2627.337 


2629.17 


-1.83 


12 0) 


2805.6 


2804.41 


1.19 


101) 


2906.074 


2905.53 


0.54 


040) 


2993.0 


2992.36 


0.64 


21) 


3092.481 


3092.90 


-0.42 


002) 


3201.433 


3202.20 


-0.77 


210) 


3364.57 


3364.02 


0.55 


130) 


3547.1 


3545.96 


1.14 


1 1 1) 


3637.843 


3635.84 


2.00 


5 0) 


3738.6 


3737.55 


1.05 


31) 


3829.34 


3829.52 


-0.18 


300) 


3922.61 


3923.95 


-1.34 


012) 


3929.12 


3928.97 


0.15 


2 2 0) 


4100.58 


4099.75 


0.83 


2 01) 


4179.938 


4177.90 


2.04 


14 0) 


4286.82 


4286.42 


0.40 


121) 


4369.1 


4366.68 


2.42 


10 2) 


4461.07 


4458.63 


2.44 


060) 


4482.57 


4481.44 


1.13 


041) 


4564.22 


4565.11 


-0.89 


310) 


4652.0 


4651.01 


0.99 


2 2) 


4656.34 


4656.38 


-0.04 


003) 


4754.209 


4753.80 


0.41 


230) 


4835.05 


4835.05 


0.00 


211) 


4905.52 


4901.97 


3.55 


15 0) 


5025.2 


5025.31 


-0.11 


131) 


5098.0 


5097.07 


0.93 


112) 


5180.54 


5177.47 


3.07 


400) 


5205.81 


5207.49 


-1.68 


070) 


5224.55 


5224.03 


0.52 


51) 


5298.16 


5299.39 


-1.23 


032) 


5377.91 


5377.30 


0.61 


3 2 0) 


5384.41 


5384.69 


-0.28 


301) 


5437.54 


5434.25 


3.29 


013) 


5469.66 


5468.64 


1.02 


2 4 0) 


5568.41 


5568.93 


-0.52 


2 21) 


5630.36 


5626.62 


3.74 


2 2) 


5701.41 


5698.37 


3.04 


16 0) 


5762.23 


5762.25 


-0.02 


141) 


5826.29 


5826.29 


0.00 


12 2) 


5898.94 


5896.50 


2.44 


410) 


5930.66 


5930.69 


-0.03 


080) 


5965.61 


5965.30 


0.31 


10 3) 


5984.705 


5981.83 


2.88 
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TABLE E. Continued 



(vi V 2 V3) 


^obs 




rriobs TTical 
£j v — £j v 


(0 61) 


6030.71 


6032.06 


-1.35 


(0 4 2) 


6101.80 


6101.79 


0.01 


(3 3 0) 


6112.11 


6112.77 


-0.66 


(311) 


6156.25 


6152.20 


4.05 


(0 2 3) 


6183.61 


6182.94 


0.67 


(0 04) 


6275.98 


6274.83 


1.15 


(2 5 0) 


6299.70 


6300.57 


-0.87 


(2 31) 


6351.40 


6350.62 


0.78 


(212) 


6414.16 


6411.21 


2.95 


(5 0) 


6475.05 


6475.27 


-0.22 


(170) 


6497.60 


6496.88 


0.72 


(151) 


6552.84 


6553.69 


-0.85 


(132) 


6616.53 


6614.85 


1.68 


(4 2 0) 


6653.54 


6654.00 


-0.46 


(401) 


6676.86 


6672.05 


4.81 


(113) 


6693.12 


6692.31 


0.81 


(0 9 0) 


6705.23 


6705.18 


0.05 


(0 71) 


6771.44 


6762.82 


8.62 


(0 5 2) 


6823.80 


6824.39 


-0.59 


(3 4 0) 


6837.75 


6838.95 


-1.2 


(3 21) 


6872.10 


6870.35 


1.75 


(0 3 3) 


6897.37 


6896.11 


1.26 


(3 2) 


6921.67 


6920.44 


1.23 


(014) 


6979.21 


6977.98 


1.23 


(2 6 0) 


7029.48 


7029.18 


0.30 


(2 41) 


7072.23 


7072.92 


-0.69 


(2 2 2) 


7125.60 


7123.90 


1.70 


(2 3) 


7192.29 


7191.82 


0.47 


(510) 


7193.35 


7192.88 


0.47 


(180) 


7231.06 


7228.80 


2.26 


(161) 


7277.83 


7278.63 


-0.80 


(142) 


7332.45 


7331.72 


0.73 


(43 0) 


7374.57 


7375.49 


-0.92 


(411) 


7386.33 


7382.08 


4.25 


(12 3) 


7403.04 


7402.08 


0.96 


(010 0) 


7443.09 


7443.40 


-0.31 


(104) 


7478.02 


7476.38 


1.64 


(0 81) 


7492.23 


7491.24 


0.99 


(0 6 2) 


7544.62 


7544.93 


-0.31 


(3 5 0) 


7562.47 


7561.88 


0.59 


(3 31) 


7587.04 


7587.25 


-0.21 


(0 4 3) 


7609.57 


7607.70 


1.87 


(312) 


7627.14 


7627.15 


-0.01 


(0 2 4) 


7681.49 


7679.79 


1.70 


(6 0) 


7730.08 


7727.56 


2.52 


(2 7 0) 


7757.29 


7753.98 


3.31 


(0 5) 


7766.28 


7766.57 


-0.29 


(2 5 1) 


7791.18 


7792.56 


— 1.38 


(2 3 2) 


7834.97 


7835.43 


-0.46 


(213) 


7888.16 


7884.58 


3.58 


(5 2 0) 


7909.46 


7905.79 


3.67 


(5 01) 


7903.54 


7909.62 


-6.08 


(19 0) 


7962.27 


7957.47 


4.80 
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TABLE E. Continued 



(Vl V2 V3) 


^obs 




TTiobs TTical 
£L/ V — Ej v 


(171) 


8000.93 


8000.40 


0.53 


(15 2) 


8046.44 


8046.39 


0.05 


(44 0) 


8093.61 


8091.13 


2.48 


(4 21) 


8093.10 


8093.50 


-0.4 


(133) 


8110.13 


8110.59 


-0.46 


(40 2) 


8120.70 


8123.66 


-2.96 


(114) 


8174.27 


8173.90 


0.37 


(0110) 


8178.27 


8179.34 


-1.07 


(0 91) 


8218.84 


8216.63 


2.21 


(0 72) 


8264.28 


8263.22 


1.06 


(3 6 0) 


8284.17 


8280.23 


3.94 


(3 41) 


8299.45 


8301.69 


-2.24 


(0 5 3) 


8320.00 


8317.42 


2.58 


(3 2 2) 


8330.35 


8333.20 


-2.85 


(3 3) 


8374.58 


8379.91 


-5.33 


(0 34) 


8382.64 


8382.98 


-0.34 


(610) 


8441.44 


8439.54 


1.90 


(015) 


8457.15 


8457.35 


-0.20 


(2 8 0) 


8482.12 


8473.96 


8.16 


(2 61) 


8507.33 


8508.48 


-1.15 


(2 4 2) 


8542.25 


8544.87 


-2.62 


(2 2 3) 


8585.54 


8584.47 


1.07 


(511) 


8608.92 


8613.47 


-4.55 


(5 3 0) 


8623.34 


8623.21 


0.13 


(2 4) 


8652.27 


8658.91 


-6.64 


(1100) 


8690.72 


8682.00 


8.72 


(181) 


8721.11 


8718.07 


3.04 


(162) 


8758.28 


8758.18 


0.10 


(4 31) 


8797.95 


8797.80 


0.15 


(4 5 0) 


8809.81 


8806.33 


3.48 


(412) 


8817.61 


8817.08 


0.53 


(143) 


8816.65 


8823.37 


-6.72 


(1 2 4) 


8868.35 


8869.19 


—0.84 


(012 0) 


8911.29 


8911.52 


-0.23 


(105) 


8941.28 


8937.70 


3.58 


(0101) 


8944.50 


8943.12 


1.38 


(70 0) 


8968.55 


8962.76 


5.79 


(0 8 2) 


8982.08 


8979.28 


2.80 
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TABLE Fa The observed vibrational band origins [cm ] of A symmetry of H2O are compared to those computed using the 
present method and to those computed in the Reference work by Lemus et al. [23] . A few predicted levels are also listed. 



(Vl V2 V3 ) 




peal 


TT'obs Trical 


EiRef 


rpobs TTiRef 


(U 1 UJ 


ioy4. to 


1 Cfi/I p. 

1594.0 


U.2 


1091.00 




6.1 


(U 2 UJ 


O 1 K 1 CO 
olol.OO 


01 ri 7 

3151. i 


n 1 
— U.l 


1 CO /to 

3152.42 


A O 

—U.o 


(1 U UJ 


ODD ( .UO 


ODOD.D 


U.O 


Qfi CQ 1 Q 

00O0.I0 


1 1 

— 1.1 


(U UJ 


4000.O 


4DD( .O 


—u.o 


ARR(\ CO 

4ooy.oz 


1 
— Z. ( 


(1 1 Uj 


OZ.30.U 


5235.4 


— U.4 


COO /) A C 

5234.45 


a a 
U.O 


(U 4 UJ 


£J1 O/l HQ 

0lo4.Uo 


OloO.Z 



— z.z 


oloo.Uo 


1 A 
l.U 


/ 1 o n\ 
(1 2 Uj 


Of f 0.1 


0/ (0.3 


1 

— 1.2 


0/8U.4/ 


— 0.4 


/'o n n\ 
(z U UJ 


70A1 K/1 

( zUl.o4 


1 zUl. < 


— U.z 


7000 79 
< ZUU. < O 


A Q 
U.O 


/'n n o\ 
(U U 2j 


TA A C H7 

f 440.U / 


7/1 /ICQ 

( 440.0 


1 
— U. 1 


*7AACl 7/1 

/ 44y. / 4 


A 1 

—4. ( 


(U 5 UJ 


7c: A O OA 

f o4z.oy 


7C A Ci C 


— ( .1 


7CQC fl/1 

<5o5.y4 


D.O 


^1 UJ 


0070 no 


8275.0 


— 1.0 


OZ 1 0.01 


/I Q 

—4.0 


/o 1 n\ 
(2 1 Uj 


oYoi.oy 


8 (63.U 


— 1.4 


07/;n 7 


A A 

u.y 


(U UJ 


QQAA C 


cool c 
ooyl.o 


— l.U 


CO 70 OC 

00 1 Z.05 


1 < . < 


/n 1 o^ 
(V 1 2 J 


nnnn i a 

yuuu.14 


yuuu.y 


n 
— U.8 


syyo.45 


7^ 


(1 4 Uj 




070fi /I 

y 1 zD.4 




0701 Ofi 

y 1 zi.yo 




(U 1 Uj 




1f\1 O A Q 

1U1o4.0 




1 m /i n 
1U14U.Z 




/o n\ 
(2 2 Uj 


1 noo A A 
!Uzo4.4 


lUzoo. ( 



— 2.3 


1 noon 1 
lU2oy.l 


—4. f 


(U z zj 


!Uoz4.o 


lUozo.l 


1 
l.z 


1 n coq 7 
lUOzo. 1 


A C 
U.O 


(3 U UJ 


lUoyy. / 


1UDU3.4 


—3. ( 


1 n cn7 n 

iu5y ( .y 


1 O 

1.0 


(1 U 2j 


1 nose n 


1 no£? 1 
IU0D8.I 


U.8 


1U8 (0.8 


1 A 

— i.y 


(,0 1 UJ 


1 1 on 

izioy.z 


1 Ol A A K 

12144.0 


—0.3 


1 01 ot? n 
12136. U 



6.1 


/1 1 o\ 
(1 1 2J 


1 O A AT a 

lz4U 1 .0 


1 /i no 7 
12408. ( 


1 -1 
— 1.1 


1 a n /i c 
12404.5 


1 
0.1 


^z 4 UJ 


1 ooac: 1 
lozUo.l 


lozUO.O 


O /I 

— U.4 


1 Q 1 OQ n 

loiyo.y 


11. z 


/n /i o\ 
(U4 2J 


1 O A A q a 

lo44o.U 


13400.0 


— / .8 


1 /i c7 n 
1345 (.U 


A A 

— y.u 


o\ 
(1 z zj 


loo4z.z 


1 Qfi /I Q O 

lo04o.U 


O.O 


1 Qfi /I Q /I 

lo04o.4 


— l.z 


/'o n o\ 
^z U zj 


1 OQOQ O 


1 Q Q 1 Q K 

loolo.O 


y.s 


la8zo.y 


1 A 

1.4 


1 o\ 
(1 z zj 


1 OA1 A A 

loyiu.y 


loyio.o 


— z.y 


1 QQ1 /I C 

loyi4.5 


a. 
— o.D 


(4 U UJ 


1 A OOI O 

14zzl.z 


1 /I 1 CI f\ 

i4ziy.u 




z.z 


1 /I O0 1 7 
14ZZ1. / 


A K 

—U.o 


(U 4J 


1 A COfi C 

140O0.0 


i4o3y.o 


— 3.U 


1 /I c 00 c 

14532.5 


/I A 

4.U 


(132) 


15107.0 


15106.4 


0.6 


15104.9 


2.1 


(212) 


15344.5 


15355.7 


-11.2 


15341.6 


2.9 


(410) 


15742.8 


15743.8 


-1.0 


15740.4 


2.4 


(2 2 2) 


16825.2 


16836.0 


-10.8 


16824.6 


0.6 


(3 2) 


16898.4 


16914.4 


-16.0 


16898.7 


-0.3 


(4 2 0) 


17227.7 


17229.8 


-2.1 


17231.2 


-3.5 


(5 0) 


17458.3 


17458.7 


-0.4 


17455.8 


2.5 


(10 4) 


17748.1 


17747.1 


1.0 


17746.5 


1.6 
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TABLE Fb The observed vibrational band origins [cm *] of B symmetry of H2O are compared to those computed using the 
present method and to those computed in the Reference work by Lemus et al. [LEMUS - final reference !!!]. 



Vi V 2 V3) 


zriobs 


peal 


TT'obs TTical 


TTlRcf 


TTiobs TTlRcf 


01) 


3755.93 


3757.0 


-1.1 


3762.48 


-6.6 


Oil) 


5331.27 


5331.8 


-0.5 


5329.2 


2.1 


21) 


6871.51 


6872.1 


-0.6 


6873.16 


-1.7 


101) 


7249.81 


7252.0 


-2.2 


7253.5 


-3.7 


31) 


8373.85 


8374.6 


-0.8 


8377.65 


-3.8 


1 1 1) 


8807.0 


8809.6 


-2.6 


8804.8 


2.2 


041) 


9833.58 


9835.0 


-1.4 


9832.29 


1.3 


121) 


10328.7 


10331.7 


-3.0 


10331.3 


-2.6 


2 01) 


10613.4 


10618.1 


-4.7 


10615.4 


-2.0 


003) 


11032.4 


11034.0 


-1.6 


11032.8 


-0.4 


51) 


- 


11246.6 


- 


11229.4 


— 


131) 


11813.2 


11816.0 


-2.8 


11815.2 


-2.0 


211) 


12151.3 


12157.1 


-5.8 


12149.2 


2.1 


013) 


12565.0 


12567.3 


-2.3 


12557.1 


7.9 


61) 


— 


12598.1 


— 


12563.1 


— 


141) 


13256.2 


13258.0 


-1.8 


13248.9 


7.3 


2 21) 


13652.7 


13659.1 


-6.4 


13655.7 


-3.0 


71) 


13835.4 


13837.7 


-2.3 


13829.5 


5.9 


2 3) 


14066.2 


14069.9 


-3.7 


14066.8 


-0.6 


103) 


14318.8 


14320.8 


-2.0 


14319.7 


-0.9 


151) 


14640.0 


14648.4 


-8.4 


14626.2 


13.8 


2 31) 


15119.0 


15125.0 


-6.0 


15118.9 


0.1 


311) 


15348.0 


15359.7 


-11.7 


15345.9 


2.1 


033) 


15534.7 


15539.6 


-4.9 


15543.9 


-9.2 


113) 


15832.8 


15837.5 


-4.7 


15825.2 


7.6 


3 21) 


16821.6 


16834.6 


-13.0 


16824.4 


-2.8 


203) 


16898.8 


16914.7 


-15.9 


16899.2 


-0.4 


12 3) 


17312.5 


17319.5 


-7.0 


17314.3 


-1.8 


401) 


17495.5 


17498.8 


-3.3 


17500. 


-4.5 
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Algebraic-matrix calculation of vibrational levels of triatomic molecules 
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and European Theoretical Spectroscopy Facility 
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(Dated: December 1, 2008) 

Q\ ' We introduce an accurate and efficient algebraic technique for the computation of the vibra- 

tional spectra of triatomic molecules, of both linear and bent equilibrium geometry. The full three- 
dimensional potential energy surface (PES), which can be based on entirely ab initio data, is pa- 
rameterized as a product Morse-cosine expansion, expressed in bond-angle internal coordinates, and 
5_l | includes explicit interactions among the local modes. We describe the stretching degrees of freedom 

in the framework of a Morse-type expansion on a suitable algebraic basis, which provides exact 
analytical expressions for the elements of a sparse Hamiltonian matrix. Likewise, we use a cosine 
power expansion on a spherical harmonics basis for the bending degree of freedom. The resulting 
l/") , matrix representation in the product space is very sparse and vibrational levels and eigenfunctions 

can be obtained by efficient diagonalization techniques. We apply this method to carbonyl sulfide 
OCS, hydrogen cyanide HCN, water H2O, and nitrogen dioxide NO2. When we base our calcula- 
tions on high-quality PESs tuned to the experimental data, the computed spectra are in very good 
agreement with the observed band origins. 



Oh 
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X 



PACS numbers: 33.20.Tp, 31.15.-p, 31.50.Bc, 03.65.Fd, 33.15.Mt 
Keywords: vibrational states, three atomic molecules, algebraic methods 



O 

u 

■ I. INTRODUCTION 

H ■ 

The present-day fast development of new spectroscopic instruments and methods allows us to measure vibrational 
states with high accuracy even in the energy region near the molecular dissociation [3, H, Q . Characterizing these 
new experimental data is especially important for instance, for understanding the dynamics of chemical reactions 
and for analyzing spectra from remote regions (e.g., upper atmosphere and/or interstellar matter), where strong 
radiation along with low pressure can stabilize molecules in extraordinarily excited states, sometimes promoting 
1 unusual reactions. All energy levels depend sensitively on the detailed shape of the molecular potential energy surface 
^\ . (PES). Their qualitative and also quantitative description calls for the development of more and more accurate 
' theoretical approaches 0, H, H, H, 0, @] • 

, Present-day status of the art program codes for the accurate calculation of vibrational spectra of small molecules 
have been developed in J. Tennyson's (TRIATOM [4], DVR3R [H] ...) and P. Jensen's (MORBID [6], TROVE |0]) 
Q\ \ groups. They are based on different approaches and find best use for different kinds of small (in particular triatomic) 
molecules. TRIATOM focuses mainly on van der Waals complexes, and employs Legendre and Laguerre polynomials 
as basis functions. The matrix elements are integrated numerically using a discrete variable representation (DVR) 
based on Gauss- Jacobi and Gauss-Laguerre quadrature for all 3 internal coordinates. The resulting spectra depend 
on the DVR grid in a non-variational fashion. The MORBID code is useful mainly for standard rigid molecules. It 
takes an advantage of Morse oscillator basis functions, for which the matrix elements of pure stretching motion are 
know analytically. However, the bending basis functions are evaluated in the framework of a numerical integration 
technique, which can be a time-consuming step. Modern fully variational approaches use a finite basis representation 
of the vibrational space: as the TROVE implementation in principle allows to perform variational calculations for general 
polyatomic molecules of arbitrary structure. However in that approach, the kinetic energy operator is represented 
only approximately as a power expansion in terms of internal coordinates. That approximation produces adequately 
accurate energy levels only rather low in energy, where the wavefunction has a limited spread around the molecular 
equilibrium position, where the power expansion is accurate. 

The present work introduces a virtually exact variational method, which joins the main advantages of (i) simple 
algebraic forms of matrix elements, (ii) completeness of the basis set, and (iii) sparseness of the resulting Hamiltonian 
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matrix. In concrete, we describe the stretching modes with a formalism based on the Morse oscillator [9(, with 
some similarity to the algebraic approaches of Ref. [Io| . The main advantage of an potential expansion in terms 
of Morse coordinates is that one can use a quantum-mechanical basis on which the matrix representation of the 
Hamiltonian operator is sparse and can be computed analytically using algebraic techniques based on generalized step 
operators. For the bending mode we choose a cosine-power expansion of the potential and a quantum mechanical 
basis of spherical-harmonic functions. We then formulate and solve numerically the multi-dimensional problem of the 
vibrations of a polyatomic molecule in a product space of the different degrees of freedom. The sparse quality of the 
individual matrices representing the one-dimensional vibrations and the interaction terms carries forward to the global 
matrix representation of the total Hamiltonian on the product Hilbert space. The eigenvalues and eigenvectors of the 
resulting sparse matrix can then be obtained numerically quite efficiently by modern exact diagonalization tools. In 
principle, this method can be applied to molecules composed by any number of atoms, even though in practice the 
size of the product Hilbert space grows exponentially in the number of dimensions, which makes this method, like all 
virtually exact methods rather unpractical for molecules with 6 or more atoms. In this study we focus on triatomics, 
but we construct the theory in terms of basic building blocks usable for successive extensions to the multi-dimensional 
potential surfaces of polyatomic molecules, which will be the subject of future investigation. 

The present approach bears some resemblance to Lie algebraic methods based on a Heisenberg formulation of 
quantum chemistry (the second quantization of the Schrodinger equation) . Several previous studies [ToL [ill IT3 , [Hj] 
use an algebraic approach where the full Hamiltonian operator is expanded in powers of C, a Casimir operator of a 
suitable Lie algebra, for example as follows: H — Eq + 2^=2 ^MC(0(2))]\ The trouble with that kind of expansion 
is that it maps to an explicit first-quantized Hamiltonian form involving an intricate mixed potential and kinetic 
contributions, including high powers of the momentum operator and unphysical products of the momentum and 
position operators. To avoid such unwanted properties we express separately the potential V and kinetic T energy 
operators. This can be done within a convenient scheme of ladder operators defined by means of a suitable factorization 
method @,[l||. 

Section [IT] sketches the model used, including the Hamiltonian, basis functions for stretching and bending modes, 
and formulas for the matrix elements. Section Hill reports and discusses the vibrational spectra of OCS, HCN, NO2, 
and H2O as obtained with the present method. Section IIVI discusses the approximations involved in the present 
method and its future extensions. 



II. THEORY 

We perform all calculations in the bond-length-angle internal coordinates, as illustrated in Fig. [1] This is the 
most suitable choice of variables for the potential expansion involving Morse functions. R\ and R2 represent the 
bond lengths between the central atom and the two end atoms, is the bending angle at the central atom. Given 
the masses m\ and m 2 of the end atoms and 7713 of the central atom, we define the two diatom reduced masses 
[i\ = m\m^l (rn\ + 7713) , \i2 = ^2^3/ (ni2 + mz). The vibrational Hamiltonian, composed by a kinetic and a potential 
part, can be expressed in atomic units (h = m e = q e = 1) as follows: 

H(Rx,R2, 0) = T(R X , R 2 ,6) + V{R U R 2 , 9) . (1) 

The standard expression for the pure vibrational kinetic energy operator in internal coordinates [l5l [l6| is 

Id 2 Id 2 



T(Ri, R2, 



2/zi dRj 2ti 2 dR\ 



L 2 

I / :i I 2cos(9 \ f d 2 n d 

cote— 



2 \fxiRl H2R2 m % R 1 R 2 J \d0 2 89 / 

sin 9 d cos 9 {Id 1 d \ cos 9 d 2 



m 3 i?ii? 2 99 m 3 \Ri 8R 2 R2 9R\ J 7113 8R\8R2 
sin# d ( 1 d 1 d \ cos 9 



(2) 



m 3 d9\R 1 dR 2 R 2 dR% J m 3 i?ii? 2 
1/1 1 2 \ 1 



H 



8 \fXiRl \i2R\ m 3 i?ii?2/ cos 2 (6>/2) 

This expression is singular at four boundary regions: at Ri — * 0, and at the bending extremes 9 = and 9 = n. 
While the former constitute no serious practical problem, since the vibrational motion avoids these regions due to the 
strongly repulsive nature of the potential energy function V there, the angular singular points are easily reached by 
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the motion: in particular 8 = n is the angular minimum energy direction in the case of linear molecules, and this may 
require explicit care to obtain a numerically stable algorithm. Specifically, the IT Z operator in the final row, defined in 
Ref. [IB], accounts for the rotations of the molecule around an axis attached to it and coincident with the molecular 
axis when the molecules reaches its linear configuration 8 — 180°: to our purpose Yi z can be replaced everywhere with 
its eigenvalue m. The final term oc IY} describes the kinetic-energy contribution of a rotational degree of freedom of 
bent molecules (and in this case we omit it), but it is needed to account for the fourth vibrational degree of freedom 
describing the m = ±1 "IT", m = ±2 "A", etc.... axially- rotating vibrational excitations of those molecules such as 
OCS and HCN which are linear in their equilibrium geometry. For simplicity, in the main text we will stick mostly 
to m = 0, while the Appendix shall deal with the general algebra describing an arbitrary integer value of m. 
We use the following parameterization of the potential energy surface: 

k 1 +k 2 +k 3 <N c 

V(Rx,R 2 ,8)= J2 a klk2k3 v 1 {R 1 pv 2 {R 2 ) k -u{e) k ^ . (3) 

fci ,k 2 ,k 3 =0 

This power expansion is realized in terms of Morse-related functions 

ViiRi) = e" a *(^- R --) - 1 , (4) 



for the stretching degrees of freedom i = 1 or 2, and in terms of trigonometric expressions u(8) — [cos 8 — cos Ami J for 
the bending degree of freedom. Similar parameterizations are employed e.g., by P. Jensen in the MORBID code 17j . 
also D. Xie and coworkers use this Morse-cosine expansion for several molecules, with the parameters first fixed to fit 
an ab-initio PES, and later adjusted to reproduce spectroscopic experimental data (l8l. Il9||. 

As customary, to compute a variational solution of the eigenproblem corresponding to H, we resort to the expansion 
of the eigenfunctions of H on a product basis: 

* v {Ri,R 2 ,8)= hh *j /^.W;. (5) 

jl ,32 ,33 

where v is a complete set of vibrational quantum numbers characterizing an eigenstate of H , cj, 1 n 33 are yet-to-be- 
dctermined expansion coefficients, and 

<S> nnn {R 1 ,R 2 ,8)=4> ln {R 1 )fo n {R 2 )Y J M- (6) 

The 1-dimensional basis functions (fiij^Ri) for the stretching modes were discussed in detail in Refs. [tl IT3 |. In 
particular, we use here the generalized quasi-number state basis [3 HO], defined as follows: 



= \ L,o lJ t m V? ^T" 1 ^) ' 3i = 0, 1, 2, . . . , i = 1 or 2 , (7) 

w 1 (2<7i + Ji) 



with 



V, 



= yi {R t ) = (2 Sl + l)e~ a ^- R —K (8) 



Here L£ are generalized Laguerre polynomials, Y is the standard Gamma function generalization of the factorial |2ll | 
and Si and Ci are suitable positive parameters, whose value we discuss below. 

For the angular variable we use the following basis derived from the spherical harmonics 



Yi{9) - s/2izY lQ {0,<p) = NtPticose) , with N t = \ l+ - (9) 

[Pi(cos8) are standard Legendre polynomials Pi(z) = -£j(z 2 — 1)V(2 Z '0]i w hich form a convenient orthonormal basis 
set over the range < 8 < n, in the measure dcosd. This basis was used successfully by Carter and Handy [TB|. In 
that work the authors evaluated the matrix-element integrals numerically by Gauss-Legendre quadrature, while here 
we use analytical expressions detailed below. 

As the Hamiltonian is a sum of terms, its matrix elements are also expressed as a sum of individual terms, and 
each one of them is expressed as a product of operators acting on the Ri, R 2 and 8 variables. Accordingly, the matrix 
elements of H are computed as sums of products of terms, each of which refers to one oscillator individually. For 
example, the matrix elements of one of the kinetic contributions in Eq. ([2]) are evaluated as 

(*AAial - ^rJ^UUO = ~^\^J0 (fel^l^> {Y j3 \co S 8\Y f3 ) . (10) 
Each of these 1-dimensional matrix elements is then evaluated using the algebraic methods described below. 
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A. Stretching-coordinate matrix elements 



As a first step, all independent terms in the Hamiltonian of Eq. |T]) must be expressed as functions of the Morse 
variable yi(Ri) of Eq. ([5]). In other words, the Vi{Ri) ki potential terms in Eq. §5§ and the R~ l and R~ 2 terms in 
the kinetic energy operator Eq. ^ must be expressed in terms of the corresponding y.^Ri). This is accomplished 
immediately for the potential terms [Tij , by choosing for each mode i the same value of cti in the potential expansion 
terms (U]) and in the definition of yi(Ri), Eq. ([5]), so that Vi = yi/(2si + 1) — 1. The matrix elements of R^ and R~ 2 
on the basis of Eq. ([7]) could be computed exactly by numerical integration, but that approach would be contrary to 
the general spirit of the present method and very inefficient, since that procedure would produce a non-sparse matrix. 
We prefer to fit the kinetic terms with a sum of powers of vf. 

Bp t 

R7P ^Y^b pi3Vl {R t y, p=l,2. (11) 

3 

An example of the quality of such a fit is illustrated in Fig. [21 The fit targets a reasonably wide region around the 
equilibrium point Ri mm , where it is extremely accurate, but it deteriorates especially in the large- Ri dissociation 
region. The fitting function could be forced to reach the same large- Ri limit as R~ 2 , but then the fit would deviate 
much more in the most important region near the minimum. Accordingly, we prefer to accept this deviation near 
dissociation, which causes negligible numerical error to the final spectra. The fitted i?i-ranges and the best-fit 
coefficients b in and buj of R~ x and R~ 2 for OCS, HCN, NO2 and H 2 are reported in the Table A of the supporting 



information [22j. Root mean square (RMS) deviations ranging from 3 to 20 cm are obtained for the fits based on 
Eq. (fTTj) . of the type exemplified in Fig.O In practice these deviations affect the computed vibrational band origins 
at a completely negligible level. 

The full details of the formalism to construct the analytical matrix elements of the functions Vi(Ri) of Eq. ([4]) are 
reported elsewhere [§, [23| . Here we just note that the idea underlying the algebraic approach is ultimately related 
to supersymmetry, which provides several analytical relations for a class of exactly-solvable problems, including the 
Schrodingcr equation for the Morse potential. We only report the basic expression of the matrix elements of the 
exponential function and of the first derivative term: 

(ikjle-"* 1 *'-*"**^') = — 5 "'^ + 2( 2 8 + /i — = 6 "' +1 > ( 12 ) 



W*' = y 



i\-^~\4>i3') - ~7T i^ij' ~ ^Jii'+l) J ( 13 ) 



where Cjj = + 2<Ji — 1). Based on these expressions, all matrix elements of every stretching term in the 

potential expansion Eq. and in the kinetic terms expressed as in Eq. (jlip can be computed exactly. The matrix 
representation of Vi(Ri) ki is (2fc; + l)-band diagonal. 

The computation of the spectrum of one pure-stretching mode through the exact diagonalization of algebraic ma- 
trices can be made substantially more efficient by choosing a properly adapted quantum-mechanical basis, specifically 
tuned to the molecular potential [l4| ■ A better convergence of the calculation of a single oscillator improves the con- 
vergence of the full 3-dimcnsional spectral calculation for the triatomic molecule. For example, for the first stretching 
oscillator we consider 

1 d 2 Nc 

Hi(Ri) = - — ^+J2 a ^ooV 1 (R 1 ) kl . (14) 



fcl=0 



The shape of the basis wavefunctions Q is tuned by four parameters i?i m i n , Q^i, si and a\. We fix i?i m ; n to the 
relevant equilibrium bond length, and in order to preserve the matrix sparseness we set a\ to the value used in the 
definition of ui(i?i). By tuning the remaining s\ and o\ one has a sufficient basis flexibility to improve substantially 
the numerical convergence of the single-oscillator problem [14] Eq. (|14p . We optimize the parameters s, and <Ji by 
the minimization of the sum of the Nj, lowest bound-state eigenenergies of the dimcr, as defined in equation (25) of 
Ref. [ijj]- We optimize the energies of Nb = 20 stretching bound states, which cover and far exceed the spectral range 
of interest. To get sub-cm -1 accuracy, a number of basis functions N<n ~ 30 in the one-oscillator basis are usually 
sufficient. The values of the final optimized parameters Sj and Cj adopted in all calculations are collected in Table [T] 
In fact, when A^j is large enough the spectral accuracy depends only weakly on the value of the basis parameters Si 
and (Tj, so that the reported values are not particularly critical. 
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B. Bending-angle matrix elements 

On the spherical harmonics basis Eq. ([9]), it is straightforward to express the matrix elements of powers of the 
cos 9 function. The bending-angle dependence of the potential energy surface is then conveniently fitted to powers of 
(cos# — cos# min ). In this basis, the bending potential matrix is sparse and can be expressed analytically. To evaluate 
the matrix elements of powers of z = cos 8, we use the following relations (2lj : 

I 7 

(YAzlYj) = Si i+ i + I Si ,_i , (15) 

WM " V(2J + 1)(2Z + 1) ' J+ vW+lP+1) ' 



+ c oil + 1) 



(j+/ + l)V(2i+l)(2/+l) ' J (j+Z + l)V(2i + l)(2Z + l) 



and in general, recursively 



Wl ' 3> V(2j + l)(2j + 3) V " ' V(2j + l)(2j-l) ' 

Like for stretching matrix elements, the matrix representing u(9) k3 is (2&3 + l)-band diagonal. 

The derivation of the kinetic bending matrix elements is described in Appendix [X] The final result is the following 
simple tridiagonal expression: 

(H|T b8nd |^) = \ (-L + -L) j(j + 1)S U (18) 



2 \niR 1 ^2-^2 
1 



m 3 RiR 2 v /(2/ + l)(2j + l) 



(Z 3 5 ij+ i +j 3 < 5 ij _i) . 



The basis functions Eq. ([5]) are independent of the shape of the bending potential, which means that, for bent 
molecules in particular, the number of basis functions required by a given degree of convergence could become 
unreasonably large. Unfortunately, there is no way of optimizing the individual functions to any given angular 
potential, as we do for the stretching functions. One way to overcome this problem is to replace the angular functions 
© with suitably optimized linear combinations thereof, for example, obtained as low-energy eigenstates of a purely 
bending 1-dimensional problem based on Tb en d plus the aoofc 3 part of the potential expansion ([3]), as was suggested 
e.g., in Ref. [l5| and Ref.s therein. Such kind of approach can lead to significant basis size reduction, but also, 
unfortunately, to entirely nonzero angular matrix elements, eventually leading to a less sparse total matrix of H. 



III. RESULTS 



We come to illustrate the application of the proposed algebraic method to the calculation of the vibrational spectra 
of real triatomic molecules. We target mainly the purely vibrational energy levels (J = 0). As trial systems we select 
two linear molecules, OCS and HCN, and two bent ones, NO2 and H2O. They are all well studied by experimental as 
well as theoretical techniques which provide accurate spectra to compare with. In particular, for PCS, NO2 and H2O 
the literature offers realistic PESs parameterized in the form §5§ suitable for our method [H, 0, Hj] . Accordingly, 
we have taken the PES function parameters ak 1 k 2 k 3 for OCS, NO2 and H2O exactly as in the referred publications 
[25l | . For HCN, we use the potential surface provided by Ref. [la] to fit an expansion of the type ([3]). All these 
PESs are determinated starting from ab initio points, but then the function parameters ak 1 k 2 k 3 are adjusted so that 
the deviation between the calculated vibrational levels and the observed spectra is minimized. In the energy region 
considered, these PESs are accurate enough to allow one to predict even highly excited vibrational levels with great 
precision. 

The values of the Morse exponential parameters aj and equilibrium positions Ri m - m are collected in Table |TTJ We 
fit the algebraic potential parameters a,k 1 k 2 k 3 for HCN to 737 points obtained using the PES function as given in 
Ref. [Ill, restricted to the HCN side of the HCN-HNC isomerization transition. The angular range is 9 = 180° to 90°, 
the coordinate ranges are J?hc / Qq = 1.4 to 4.4 and Rcn/ao = 1.6 to 3.0, with energy up to 45 000cm~ 1 . Table B in 
the supporting information [22| reports the resulting best fit parameters Ofc 1 fc 2 fc 3 , which produce a RMS deviation of 
29 cm- 1 . 
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As the bending frequency is usually smaller than stretching, and as the bending basis ([9]) cannot be optimized to 
the problem at hand, it is generally necessary to include a larger number of bending basis states. The minimum 
size iVd = iVd 1 A'd 2 3 of the matrix for the total Hamiltonian ([I]) which needs to be diagonalized for a very accurate 
and reasonably stable spectrum is eventually very moderate, of order Na ~ 10 4 . The nonzero matrix elements of H 
are of the order of ~ [| max(iV c , B\ 1, B\ 2, B2 1, B22) + l] 3 m each row of the global matrix. For Bki = 5 > N c , as in 
most calculation of the present work, this yields approximately 80 nonzero matrix elements per row, i.e. a very sparse 
matrix to store and diagonalize. 

Once the full matrix is constructed and stored, it is diagonalized using the Jacobi-Davidson method provided by 
the PRIMME package [27ll28| for sparse matrices. Table Hill summarizes the quality of the computed spectra, in terms of 
root mean squared (RMS) deviations between the observed and calculated vibrational energy levels. These levels cover 
an energy range extending from the ground state up to several thousand wavenumbers, as indicated in Table IIIII The 
good agreement with experiment shown by the small RMS deviations indicates both the high quality of the considered 
PES parameterizations and the satisfactory accuracy of the employed method. The complete spectral levels obtained 
by means of the present approach and their detailed comparison to experimental data for OCS, HCN, NO2 and H2O 
are collected in Tables C, D, E, F in the supporting information [22j |. For OCS we compute and compare to experiment 
also numerous II (m — ±1) states. 

For the molecules for which spectra were computed before based on the same PES (but with a different approach 
to the solution of the quantum-mechanical problem) we obtain an essentially equivalent accuracy of the spectra, as 
shown by the similar RMS deviations reported in Table IIIII Basically all discrepancies with exp eriment are to be 
attributed to the lack of accuracy of the adopted PES. Table D in the supporting information [22j collects also a few 
predicted energy levels of HCN in the 7000 to 12000 cm -1 spectral range. These calculated values are also in good 
agreement with previous calculations by Mourik et al. Table F in the supporting information [22j reports also a few 
predicted levels for H 2 0. 

The level assignments in terms of local vibrations is not always straightforward. The wave function of each vibra- 
tional state can be analyzed and it always results in a complicated admixture of excitations of all three local vibrational 
modes. For instance, the OCS vibrational excitation at E v = 2937.2 cm -1 consists mainly (35%) of v = (140) plus 
21% v = (0 6 0), plus other minor components, which confirms the traditional assignment (14 0) reported in Table C. 
However, for example the largest local component, v = (0 8 0), of the OCS vibrational excitation near 3990 cm -1 is 
26%, while the component on its standard local assignment v = (16 0) is smaller (18%). 

IV. DISCUSSION AND CONCLUSIONS 

In the present work, we demonstrate the possibility to compute the vibrational spectrum of an arbitrary triatomic 
molecule based on its PES and using algebraic techniques for the analytical determination of the matrix elements. 
The advantages of this method include simple formulas for the matrix elements, moderate total Hilbert space size, 
sparseness of the Hamiltonian matrix, all cooperating to a fast and efficient diagonalization. The generalization of this 
method to four-atomic and larger molecules is in principle straightforward, and in these higher-dimensional contexts 
the reduced basis size for each degree of freedom is even more crucial. These advantages make this method a very 
promising tool for the analysis of the vibrational levels of small molecules. 

The input of the present approach includes only the atomic masses and a reasonably dense numerical sampling of 
the PES: this is in principle within reach of ab-initio electronic structure quantum chemical calculations. We did carry 
out a calculation of the spectrum of HCN based on the ab-initio PES generated by a large set of DFT-LDA (density 
functional theory in the local-density approximation) calculations. Due to the known drawbacks of the LDA, the 
resulting spectrum shows deviations from the experimental levels of a few hundred cm -1 , but it proves the possibility 
to compute even highly excited vibrational molecular states entirely from first principles, with no parameters adjusted 
to the observed spectroscopical data. When highly accurate quantum chemical methods are applied to a relatively fine 
and extensive determination of the PES of a triatomic molecule, it will be possible to obtain much better predictive 
power for a calculation entirely free from any experimental input. 

If the PES is know to a high degree of accuracy even in the energy region near and slightly above dissociation, 
the present method can take advantage of the possibility to extend the basis set beyond the number of bound 
states to describe reliably weakly-bound pre-dissociation states in the quasi continuum. In the future, the present 
algebraic approach could be extended to the study of resonances in the continuum, bound-to- free transitions in infrared 
absorption, Franck-Condon processes and in principle even atom-molecule and molecule-molecule collisions. Other 
perspective applications include the area of non-rigid molecules (van der Waals complexes, quasi-linear molecules) 
[ll| and potentials with many minima such as those occurring in torsional oscillations [12J . Note finally that a reliable 
description of the quantum vibrational dynamics could allow us to exploit the manageable Hamiltonian structure 
produced by the present method to the study of intramolecular vibrational-energy redistribution. 
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APPENDIX A: ANGULAR MATRIX ELEMENTS 



The kinetic energy operator, Eq. ([2]), contains terms proportional to cot (9 and to [cos(#/2)]~ 2 , which are singular 
at the boundary points 6 = and tt. In fact 9 = n is an especially important configuration for linear molecules, and 
it should then be described smoothly. Luckily, since the integration implied in the matrix elements is performed over 
the variable z = cos 9, this divergence is not a problem, and all relevant matrix elements can be computed analytically 
using several properties of the spherical harmonics. For brevity, we replace the ^-independent kinetic coefficients with 
the shorthand 



a i 



a 2 = 



--( — 

2 \ixiRl 

1 

m 3 RiR 2 ' 



1 



and define the pure bending part of T as 

lbend(0) = ( a l+02C0S( 



d9 2 



cot 



89 



— a 2 sin 



89 



a 2 — 0,1 



1 



1 + cos 9 



In terms of z = cos ( 



a 



d cos 8 



(l^ 2 ) 1/2 |,and 



do 2 



COS ( 



-COS0 



(Al) 
(A2) 

(A3) 



,2\ d 2 



the angular kinetic energy is conveniently decomposed as 

7bcnd(z) = ai(l - ^ 2 )ap 
+ a 2 z(l - z 2 ) J 
+ ai (-2z&) 

ld_ 

dz / 
1 i)z 



d 

d cos 8 



<>::- 



o 2 (-2z 2 £) 

«2(1-Z 2 )f 



a 2 z 



Ti 
T 2 
T 3 
Ti 
T 5 
T 6 
T 7 . 



(A4) 



We derive the matrix elements of Tb en d in the general case of arbitrary integer m. Because for the m ^ "II" , "A" , 
etc. vibrational excitations of linear molecules, the Tj term is singular at 9 = 180°, all its matrix elements in the Yiq 
basis diverge: for general m we need then to consider the natural extension of the basis of Eq. ([9]), namely the one 
provided by the Yi m spherical harmonics: 



Yr(9) = V2^Y lm (9,Q) = N lm Pr (cos 9), 



with Ni m = 



(2l + l)(l-\m\)\ 



2(Z+|m|)! 



1/2 



(A5) 
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for I > |m|. For compactness, in the following we write the matrix elements in terms of associated Legendre polynomials 
PJ n {z) = {— l)(l m l+ m )/ 2 (1 — z 2 )l m l/ 2 dPi(z)/dz\ m \, i.e. we omit the normalization factor iV im , and we use the notation 



d pm 
dz r l 



P™. In the derivation of the Tb on d matrix elements, we make use of the following identities for the associated 



Legendre polynomials: 



(j + 1 - \m\)Pp +1 (z) = (2j + l)zP j (z)-(j + \m\)Pp_ 1 (z), 



■3+1 

(l-z 2 )Pp'(z) 



{l-z 2 )Pf'{z) 



-jzP™(z) + (j + \m\)P™(z), 



2zPf{z) 



J(J + 1) 



P™(z) . 



As a direct consequence of Eq. (|A6|) , Eq. ([T5|l generalizes to 



(Yr\z\Y" 



^ } N, 

and the recursive relation (|17[) generalizes to 

.7 + 1 



I m 
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2j + 1 



i + H 

2j + l 



(^"l^iy/ 1 ) 



2j + l 



AT 
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j + TO 



j + 1 m 
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7-lm 



These relations are also useful for the calculation of the potential matrix elements. 
The individual kinetic terms are then expressed as 



(pr\Tx\pp) = ai (pr\(i-z 2 )\pp") 
{prmpj 1 ) 



ai {PT\2z\Pf) - ad(j + i)(PT\PD + aim 2 (Pr\ T 



a 2 (Pr\z(l-z 2 )\P™"} 



(pr\T 3 \pp) 
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■\ p " 



-a 2 (Pp\2z 2 \Pp') 

a 2 (PP\(l-z 2 )\Pp')=-a 2 j(Pr\z\P? 
-a 2 (Pr\z\P; n ) 

\PP). 



a 2 (j 



i\)(pr\pr-i) 



a 2 -oi 2 

■ TO (P, 



2 " 1 '1 + z- 

By collecting everything together and simplifying we obtain 

(pr\T hend \pr) = - ai j(j + i)(priFr>-«2(j + i) 2 (pn^i^r> 



-a 2 U + H){Pr\P£i) 



ai + a 2 



(Pi" 



1 



3-1* ■ 2 - ■ 1 jg; 

Accordingly, the bending kinetic matrix elements in the normalized basis Yj m are 
(Yr\T hcnd \Y.^ ■ • " ; - " ; " J 



a 1 j{ ] + l)8 l3 -a 2 (j + lYWzlYj) 

-« 2 (j + imi)#2- ^ + to 2 (irij^iiy 



+ 



"I ./(./ + 1 ) <>/., 
02 

' V(27TT)(27TT) 

ai + a 2 2 



to 2 Sij +1 +j 2 \/j 2 ~ 



m l (YP\- \Y™) 

1 — z J 



(A6) 
(A7) 

(A8) 



(A9) 
(A10) 

(All) 
(A12) 

(A13) 
(A14) 
(A15) 
(A16) 

(A17) 
(A18) 



(A19) 



where we have inserted the matrix elements of z from Eq. (|A9[) . and combined the two non-Hermitian terms propor- 
tional to a 2 into a final explicitly Hermitian expression. The important special case to = is reported as Eq. (|18|) in 
the main text. The matrix clement in the final term is 



{Y[ | ^ j Yj ) to to ^max [^max ~t~ 1] -^-m (^max) ; 



with L 



max(l,j) . 



(A20) 
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and with 



K ±1 (l) 



1 



K± 3 (l) 



K± 2 (l) 



-(I 2 + 1-2) 

i(Z 4 + 22 3 -7Z 2 -8Z + 12). 



Unfortunately the last term in Eq. (|A19[) . the one proportional to m 2 and originated by T-j plus parts of T\ and T 2 , 
produces a matrix which, according to Eq. (IA20|) . is not sparse on this basis. This leads to significant numerical 
overhead in the calculation of the m/0 states, with respect to the m = ones. 
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TABLE I: Optimized parameters Si and at, number of target bound states Nt, and size N^i of the quantum basis for each 
stretching oscillator of all molecules studied in the present work. 



stretching 




optimized 




number of 


number of 


dimer 




parameters 




bound states 


basis functions 




s 




a 


N b 


N di 


OC 


68. 




40 


20 


30 


cs 


76. 




52 


20 


30 


HC 


30. 




0.20 


20 


40 


CN 


74. 




46 


20 


30 


OH 


26. 




0.01 


20 


30 



NO 33. 5.04 20 30 



TABLE II: The values of several relevant potential parameters for the molecules studied in the present work. The HCN 
parameters are obtained through a best fit to 737 points on the PES as computed in Ref. [26l |. 



system 


Rl min [do] 


i?2mtn [«o] 


0min [degree] 


ai [a x ] 


a 2 [a ] 


Reference 


OCS 


2.1849 (OC) 


2.9506 (CS) 


180.0 


1.2382 


1.0318 


[18] 


HCN 


2.0135 (HC) 


2.1793 (CN) 


180.0 


0.9727 


1.2290 


[26] 


N0 2 


2.2435 


2.2435 


133.767 


1.6853 


1.6853 


f24] 


H 2 


1.8112 


1.8112 


104.440 


1.1769 


1.1769 


[19] 



TABLE III: Basis size for the individual oscillators N di , total basis size Nd, and standard deviations of the computed spectra 
to the available experimental vibrational levels and to previous calculations of the same levels. 



system 


iVdi 


number of 
basis function 
N d2 N d3 


N d 


num. of 

energy 
levels 


up to 
energy 
[cm" 1 ] 


RMS deviation 
obs.-cal obs.-Ref. 
[cm" 1 ] [cm" 1 ] 


OCS 


30 


30 


50 


45xlO a 


145 


8057 


0.31 


0.26 [18] 


HCN 


40 


30 


50 


60xl0 3 


34 


12389 


14. 


12. [26] 


N0 2 


30 


30 


55 


49.5xl0 3 


143 


8979 


1.7 


2.1 [24] 


H 2 


30 


30 


50 


45xl0 3 


69 


21247 


5.6 


1.2 [19] 
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/T7-I 

FIG. 1: The scheme of internal coordinates used in the present work. 




R [a.u.] 



FIG. 2: Comparison between one of the R~ 2 function of Eq. solid line, indicated simply as R~ 2 , and its fit in terms of 
Morse-related functions, Eq. dashed line. In this calculation, relevant for the H2O stretching coordinate R, the parameters 
a ~ 1.18 fig , -Rmin — 1.81 fflo, and the expansion extends to order B2 = 5. The coefficients &2 j are adjusted to obtain a best fit 
of equally spaced points in the 1.35 to 3.2 ao range of R highlighted in the inset. 



